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ABSTRACT 

We treat the production of neutrons, photons, and neutrinos through photomeson 
interactions of relativistic protons with ambient photons in the compact inner jets of 
blazars. Internal synchrotron and external isotropic radiation due to scattered opti- 
cal/UV accretion-disk radiation are considered as target photon fields. Protons are 
assumed to be accelerated to a maximum energy limited by the size scale and magnetic 
field of the jet, and by competing energy losses. We characterize the conditions when 
the photomeson interactions of ultrarelativistic protons become effective, and show that 
the presence of the external radiation field makes possible strong energy losses already 
for protons with energies Ep > 10^^ eV. Without this component, effective energy losses 
of protons begin at Ep > 10^^ eV, and would rapidly disappear with expansion of the 
blob. 

We develop a model describing the production and escape of neutrons from a comov- 
ing spherical blob, which continue to interact with the ambient external radiation field 
on the parsec-scale broad line region (BLR). Neutrons may carry ~ 10% of the overall 
energy of the accelerated protons with Ep > 10^^ eV outside the BLR. Ultra-high energy 
gamma rays produced by photomeson interaction of neutrons outside the blob can also 
escape the BLR. The escaping neutrons, gamma rays, and neutrinos form a collimated 
neutral beam with a characteristic opening angle ^ ~ l/L, where F is the bulk Lorentz 
factor of the inner jet. Energy and momentum is deposited in the extended jet from the 
decay of neutrons at distances ld{En) ~ {En/W^''^ eV) kpc, and through pair-production 
attenuation of gamma rays with energies E^ > 10^^ eV which propagate to ~ 10-100 
kpc distances. In this scenario, neutral beams of ultra-high energy gamma rays and 
neutrons can be the reason for straight extended jets such as in Pictor A. Fluxes of 
neutrinos detectable with km-scale neutrino telescopes are predicted from flat spectrum 
radio quasars such as 3C 279. 
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1. Introduction 

Multiwavelength observations of flares from blazars, particularly in the 7-ray domain, have 
convincingly demonstrated that the compact inner jets of blazars are effective accelerators of par- 
ticles to very high energies (Hartman et al. 1999; Weekes 2000). Analyses of correlated X-ray and 
TeV gamma-ray flares in BL Lac objects lend support to leptonic models (Mastichiadis & Kirk 
1997; Catanese et al. 1997; Plan et al. 1998; Tavecchio et al. 2001; Krawczynski et al. 2002), which 
imply efficient acceleration of relativistic electrons in these sources, probably due to relativistic 
shocks. Associated acceleration of hadrons is expected with at least the same efficiency as that of 
the leptons, except perhaps for electron-positron pair jet models where few hadrons are present. 
Comparison of the radio lobe and inner jet powers indicates that jets are composed mainly of elec- 
trons and protons (Celotti Sz Fabian 1993), so that a nonthermal hadronic component is expected. 

Acceleration of hadrons in blazar jets could be directly confirmed with the detection of neutri- 
nos, provided that there are significant interactions of accelerated hadrons with ambient material 
or photon fields. Detectable synchrotron emission at TeV energies could be radiated by ultrarela- 
tivistic protons and ions, but this requires extremely strong magnetic fields ~ 20-100 G in BL Lac 
objects (Aharonian 2000; Miicke & Protheroe 2001). All other observable consequences of hadron 
acceleration result from the same interactions that produce neutrinos. These include X-ray and 
gamma-ray production from secondary leptons and gamma rays formed in hadronic interactions 
and, as we show here, deposition of energy transported by the escaping neutral radiations far from 
the inner jet. 

Models invoking interactions with ambient matter (Beall & Bednarek 1999; Pohl & Schlickeiser 
2000; Schuster, Pohl, &; Schlickeiser 2002) require mass-loaded jets. If the nuclear interaction 
energy-loss time scale t'pp = {n',^cKpp(jpp)~^ is to be less than the variability time scale t'^ari ^^^^ 
high plasma densities n'p ^> 10^/[t^ar(d)(5/10)] cm~^ are required, where n'p is the comoving thermal 
proton density, ti,ar(d) is the observed variability time scale in days, 6 is the Doppler factor, cjpp = 
30 mb is the nuclear interaction cross section, and K^p ~ 0.5 is the inelasticity. (Henceforth 
primes denote quantities in the comoving frame.) Such models would, however, be inefficient if the 
sources are to be optically thin, as is required for nonthermal X-ray escape. For Thomson-thin jets 
(tsc = npG^R' < 1), the nuclear interaction time scale 

where R' is the comoving blob radius, t'^y^ = R'/c is the dynamical (or light crossing) time scale, 
and (Jt = 665 mb is the Thomson cross sections. Thus protons can only lose < 2-3% of their energy 
on the dynamical time scale. Moreover, as shown in Appendix A, nuclear pp interaction models 
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require large masses and kinetic energies. 

A second group of hadronic models is based upon photomeson interactions of relativistic 
hadrons with ambient photon fields in the jet. Most of the models of this type take into account 
collisions of high-energy protons with the internal synchrotron photons (Mannheim and Biermann 
1992; Mannheim 1993; Miicke et al. 2002), while others also take into account external radiation 
that originates either directly from the accretion disk (Bednarek and Protheroe 1999) or from disk 
radiation that is scattered by surrounding clouds to form a quasi-isotropic radiation field (Atoyan 
& Dermer 2001). BL Lac objects have weak emission lines, so in these sources the dominant soft 
photon field is thought to be the internal synchrotron emission. Strong optical emission lines from 
the illumination of broad-line region (BLR) clouds in flat spectrum radio quasars (FSRQs) reveal 
bright accretion-disk and scattered disk radiation in the inner regions (Netzer 1990). 

In the case of internal synchrotron radiation, the energy output of secondary particles formed 
in photohadronic processes is generally peaked in the energy range from ~ lO^'^^-lO"'^^ cV in cither 
low- or high-frequency peaked BL Lac objects (Miicke et al. 2002), which implies that such models 
can only be efficient if protons are accelerated to even higher energies. This demand upon proton 
acceleration for efficient photomeson production on the internal synchrotron photons also holds 
in FSRQs, which have similar nonthermal soft radiation spectra as low-frequency peaked BL Lac 
objects. As shown by Atoyan & Dermer (2001), however, the presence of the isotropic external 
radiation field in the vicinity of the jets of FSRQs strongly improves the photomeson production 
efficiency and relaxes the very high minimum proton energies needed for efficient production of 
secondaries. The existence of a strong external radiation field is required to explain (Dermer & 
Schlickeiser 1993, 2002; Sikora et al. 1994, 2001; Dermer, Sturner, & Schlickeiser 1997; Bottcher 
2000) the luminous 100 MeV - GeV gamma-ray emission observed with the EGRET instrument on 
the Compton Gamma Ray Observatory (Hartman et al. 1996). 

In the model of Atoyan & Dermer (2001), protons are assumed to be accelerated in an out- 
flowing plasma blob moving with bulk Lorentz factor F along the symmetry axis of the accrction- 
disk/jet system. The relativistic protons are assumed to have an isotropic pitch-angle distribution 
in the comoving frame of a plasma blob, within which is entrained a tangled magnetic field. In 
our study, we determined the intensity of the internal radiation fields based on observations of 3C 
279 during the flaring state in 1996 (Wehrle et al. 1998), and calculated the high-energy neutrino 
flux expected under the assumption that the power to accelerate relativistic protons was equal to 
the power injected into nonthermal electrons which explains the observed gamma-ray emission. 
The presence of a quasi-isotropic external radiation field enhances the neutrino detection rate by 
an order-of-magnitude or more over the case where the field is absent, so we predict that FSRQs 
can be detected with km-scale neutrino detectors, whereas BL Lac objects are not as promising for 
neutrino detection. The model also takes into account the effects of relativistic neutron production, 
which can escape from the blob unless they are converted back to protons due to decay or further 
photohadronic collisions inside the blob. 
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In this paper, we further develop this model. Details of the theory are presented that describe 
the evolution of rclativistic protons, taking into account photohadronic energy losses and neutron 
escape from the relativistically moving blob. Production of neutrons has been considered earlier 
in the stationary cores of AGNs (Eichler & Wiita 1978; Sikora et al. 1989; Kirk & Mastichiadis 
1989; Bcgclman et al. 1990; Giovanoni k, Kazanas 1990; Atoyan 1992a,b), though not in the context 
of a jet model. The neutrons escaping from a relativistic blob form a beam with opening angle 
= r~^. These neutrons are subject to decay and further photohadronic interactions while passing 
through the quasi-isotropic external radiation field on the parsec scale, which we associate with 
the broad emission-line region (BLR). Ultra-high energy gamma rays produced from the decay 
of secondary pious outside the blob can escape the BLR because they are no longer subject to 
strong 77 interactions with the synchrotron photons inside the blob when they are at a distance 
of a few xR' from the blob. This results in a neutral beam of ultra-high energy (UHE) neutrons 
and gamma rays that can transport the energy to large distances from the central engine without 
significant interactions with the ambient medium until the neutrons decay or gamma rays are 
converted to electron-positron pairs through pair attenuation with diffuse radiation fields. This 
scenario naturally explains the appearance of large-scale jets that are colinear with the inner jets and 
appear straight on scales of hundreds of kpc. We also consider whether the different morphologies 
of FRI and FRII radio galaxies are a consequence of the different neutral beam properties formed 
in BL Lac objects and FSRQs. 

In Section 2, we describe the model, our treatment of photomeson production and the energy 
spectra of secondaries, and our method to estimate the radiation and magnetic fields in the blob. 
The results of the calculations are presented in Section 3, and a discussion and summary are given in 
Section 4. A comparison of neutrino production through photomeson p'j and nuclear pp interactions 
in jets is given in Appendix A, and the conditions when the direct accretion-disk radiation field is 
less important than the scattered radiation field are derived in Appendix B. 

2. Model 

Here we consider the production of secondaries formed through photomeson interactions be- 
tween relativistic protons and photons of both the internal synchrotron and the external scattered 
accretion-disk radiation fields. The scattered accretion-disk field is more important than the direct 
accretion disk radiation field when the blob is at a distance of hundreds of Schwarzschild radii 
from the central engine because, unlike the case for Compton scattering, the low-energy threshold 
for the photomeson process requires higher energy photons in the rest frame of the proton. The 
effects of directionality from photons originating from the accretion disk make this component less 
important for photomeson production than the scattered radiation field. The energy threshold is 
difficult to satisfy for photons originating from the luminous central regions of the accretion disk, 
which enter the plasma blob in a tail-on geometry. The accretion disk photons that enter from 
moderate angles originate from cooler portions of the disk at large disk radii, so that higher energy 
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protons are required to exceed the photomeson threshold energy (Bednarek and Protheroe 1999). 
A quantitative estimate of the relative importance of the direct and scattered radiation fields is 
derived in Appendix B. 

For an isotropic spatial distribution of ultrarelativistic neutrons in the frame of a blob moving 
with Lorentz factor F = \j \J\ — ]P' ^ 1, the pattern of escaping neutrons in the stationary 
frame is beamed with a characteristic opening angle Q ~ 1/F. The beaming of energetic neutrons 
would become even more pronounced when one takes into account that the angular distribution 
of external radiation is highly beamed in the direction opposite to the jet in the blob frame. For 
isotropic protons, this will give significant preference to head-on p7 collisions of the protons in the 
direction of the jet rather than in the opposite direction. Therefore, similar to the case of Compton 
scattering (Dcrmcr 1995), though possibly even more pronounced due to threshold effects, the 
external radiation field will cause a significant directional enhancement of the high-energy neutral 
beams in the jet direction already in the comoving frame. The neglect of this effect does not change, 
however, the mean energy-loss rate of protons in the blob, which are assumed to be isotropically 
distributed. It will have a noticeable impact only on the angular distribution of secondaries in 
near-threshold interactions. We have therefore neglected this effect in the paper, and a precise 
treatment of this effect is deferred to later work. 

In the following subsections, we present our method of calculations of photomeson and neutron 
production, proton energy-loss time scales and the relevant radiation fields, and spectra of secondary 
leptons, gamma rays, and neutrinos resulting from the decay of pious. Following this, we describe 
how the target photon fields are derived. We then present our treatment of proton evolution using a 
continuity-equation approach, and develop a theory for the neutron escape and subsequent neutron 
and 7-ray interactions outside the plasma blob. 

2.1. Production of Secondaries in Photomeson Collisions 

Photomeson interactions of relativistic protons with photons of the ambient radiation field 
produce secondary electrons and positrons (referred to henceforth as electrons, unless the distinction 
is necessary), gamma rays, and neutrinos primarily through the production and decay of pions. The 
relevant decay chains are 7r° — > 2 7, and tt^ ^ //^ -|- — > -h 2 -h (we do not distinguish 
between v^^ and P^). It is important that in about half of these inelastic collisions, the primary 
relativistic proton will be converted to a relativistic neutron. This is true for the single-pion 
production channel p + 7 — > n -\- tt"*" , which has approximately the same cross section as the 
channel p -|- 7 — > p -|- tt^ without neutron production, as well as for multipion production channels. 

Energy losses of relativistic protons (and neutrons) are calculated on the basis of standard 
expressions (e.g., Berezinskii & Grigoreva (1988)) for the cooling time of relativistic protons, = 
Ep/\dEp/ dt\, due to photopion production in p7 collisions. If the ambient photons have spectral 
density nph(e, fi) in the direction then the cooling time of a proton in the direction fip is equal 
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to 



(2) 







where = €jp{l — [5p cos V') is the photon energy in the rest frame of the proton, 7^, /?pC, and V' are 
the proton's Lorentz factor, speed, and colUsion angle with the photon, respectively, and K-p^{e^) 
is the inelasticity of the interaction. The threshold photon energy for pion production in and 
717 reactions in the rest frame of the incident nucleon is eth ~ 150 MeV. This means that protons 
with, say, 1 PeV energy will effectively interact with soft X-ray photons with energy > 0.15 keV. 

In the case of an isotropic photon distribution where nph(e, = n(e)/47r, equation (2) reduces 
to the form 



(Stecker 1968; Berezinskii et al. 1990). The P7 collision rate Vp^{Ep) is given by a similar expression 
as equations (2) or (3), but with the inelasticity coefficient Kp^ being dropped from the integral. 
The mean inelasticity Kp^(jp) = 1 /\tp^{'-^p)vpy{Ep)\. 

A detailed recent study of this photohadronic process is given by Miicke et al. (1999). To 
simplify calculations, we approximate cr(er) as a sum of 2 step-functions cri(er) and 0"2(er) for the 
singlc-pion and multi-pion production channels, respectively, with ai = 340 /ib for 200 MeV < Cr < 
500 MeV and cJi = outside this region, whereas C72 = 120 /xb at > 500 McV. The inelasticity 
in the single pion channel is approximated as K-p^ = K\ ~ 0.2, whereas Kp^ = K2 ~ 0.6 for 
energies above 500 MeV. Our treatment gives energy-loss pathlengths of 154, 22.0, and 14.1 Mpc 
for the photomeson interactions of ultra-relativistic protons at energies of 10^°, 3x10^^, and 10^^ 
eV, respectively, with the cosmic microwave background radiation field. These values are within 
2% of the values calculated by Stanev et al. (2000). Even at 7x10^*^ cV, where the interactions 
occur with the exponential Wien tale of the blackbody spectrum, our treatment gives an energy-loss 
pathlength of 650 Mpc, which is within 10% of the more accurate values calculated by Stanev et 
al. (2000). This approach also works well, within a few percent, for a broad power-law distribution 
of field photons u'^^i^e') oc (e')~" for different spectral indices a, and readily explains the significant 
increase in the mean inelasticity of incident protons (or neutrons) from Kp^ ~ 0.2 for steep photon 
spectra with a > 1, to Kp^ — > 0.6 for hard spectra with a < 1, giving results in good agreement 
with the treatment of Miicke et al. (1999). 

The spectra of secondary Tr'^'^-decay particles {v, 7, e) are calculated in the (5-function ap- 
proximation for the energies of the secondaries. To correctly apply the (5-function approximation, 
one has to properly take into account the different inelasticities of the single-pion and multi-pion 
production channels. In the single-pion production channel, the probability for the conversion of 
a proton to a neutron with the emission of a 7r+-meson is given by = 0.5. Given that the 
mean energies En and Ee of secondary electrons and neutrinos, respectively, are approximately 
the same, and that the vr"*" takes on average KiEp of the initial proton energy Ep^ we have that 
Ey = Ee = 0.05£^p. In the single-pion channel, a 7r°-meson is produced with the probability 1 — Cpm 




(3) 
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which results in two gamma-rays each with mean energy Ej = O.lEp after the 7r° decay. 

In the multi-pion channel, we treat the coefficient of the inelasticity K2 ~ 0.6 as a result of the 
production of three leading high-energy pions, equally distributed between tt"*", tt", and 7r° mesons, 
which becomes increasingly valid in the limit of large multiplicities. These three leading pions are 
assumed to take most of the energy lost by the primary proton, so that the remaining pions take 
away a negligible fraction of the initial energy. In this case, the mean energy carried by each pion 
is equal to 0.2Ep, and the mean energies carried by the secondary electrons, neutrinos, and gamma 
rays are the same as in the single-pion channels. The production spectra of the secondaries can 
now be found in the (5-function approximation for the differential cross sections for parent pion 
production, given the mean numbers of pions produced per pj interaction. 

The mean inelasticity Kjrf = PiKi + {1 — pi)K2 for photomeson interactions of a proton with 
energy Ep depends on the interaction probabilities pi and p2 = 1 — pi of the proton with the given 
photon field via single-pion or multi-pion channels, respectively. This relation can be inverted to 
derive the value of pi = pi{Ep), given by 

Pi = ^^^^^. (4) 

Using the numbers of neutral and charged pions produced in single-pion and multi-pion chan- 
nels, we can deduce the mean numbers of high-energy pions produced per p7 collision at proton 
energy Ep. For production we thus have n(7r°) = pi (1 — ^pn) +P2 = 1 —Pi^pn, and for the mean 
total number of charged pions n(7r+ -|-7r~) = 2p2 +PiCpn- Taking further into account that 2 7-rays 
are produced in 7r° decay, and that 3 neutrinos and 1 electron are produced in the decay of every 
charged pion, we arrive at the following formulae for the production spectra Qx{E) = dNx{E)/dt 
of gamma rays, electrons and neutrinos, given by 

Q^{E^) ^ 20 (1 - ^pnPi) Up^ilOE^) Np{10E^) (5) 

Qe{Ee) = 20{pi^pn + 2p2) iypji20Ee)Np{20Ee) (6) 
The production spectrum of neutrinos Qi,{Ei,) = 3Qe{Ee), with E^ = E^. 



2.2. Photon and Magnetic Fields 

The spectral energy density of the synchrotron radiation field ^^(e'), which is assumed to be 
isotropic in the comoving frame, depends on knowledge of the size R' of the emitting region and 
the uFi, spectral fluxes /s(e) (erg cm~^ s~^) measured by an observer during a blazar flare. We use 
the variability time scale to estimate R' ~ cSt^ar/i^ + z) although, strictly speaking, this relation 
gives only an upper limit to the source size. The comoving synchrotron photon energy density is 
then found through the relation 
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where e = Se'/{1 + z) relates photon energies in the observer and comoving frames, and is 
the luminosity distance (see, e.g., Totani (1999)). In our calculations, we consider an = 0.3, 
r^A = 0.7 cosmology with a Hubble constant of 65 km s~^ Mpc~^. 

For the spectrum of the external UV radiation field, we use a Shakura and Sunyaev (1973) 
optically-thick accretion disk radiation field that is scattered by BLR clouds. The specific specral 
energy density of the direct disk radiation along the jet axis in this model is given by (Dermer & 
Schlickeiser 2002) 

where Lad is the total luminosity of the accretion disk, h = H/Ry is the height above the accretion 
disk in units of gravitational radius Rg = GM^h/c^ = 1.5 x 10^^ M^h/MQ cm for a central black hole 
with mass M^/j. This expression employs a monochromatic approximation for the mean photon 
energy e = e*/r^/^, where e* ^ 100 eV depends weakly on Mi,h and accretion rate (Dermer & 
Schlickeiser 1993, 2002), and r = hianijj is the accretion disk radius (in units of Rg) wherefrom 
originate the photons e crossing the jet axis at an angle tp. The maximum photon energy emax is 
given in this approximation by £*/r- , where is the innermost radius of the blackbody accretion 
disk, and > 6 for a Schwarzschild metric. 

The spectrum of the isotropic component depends on the integrated disk field, which has 
spectral luminosity Ladi^) oc e^^"^ cxp{—e/emax) from the far IR up to emax- The quasar 3C 273 
displays a pronounced UV bump peaking above « 10 eV (Lichti et al. 1995). A thermal emission 
component is detected in a low state from the blazar 3C 279 at a temperature of ~ 20000 K, corre- 
sponding to emax ~ 5 eV (Plan et al. 1999). This value, however, may underestimate the effective 
maximum temperature due to uncertainties in the modeling and subtraction of the nonthermal 
component. Moreover, in a flaring state, 3C 279 may have a higher accretion rate and effective 
accretion-disk temperature. For calculations of the spectra of a typical FR II, we take emax = 20 
eV. For ~ 100 eV this implies an effective disk innermost radius ~ 8. 

The energy density of the quasi-isotropic component of the surrounding radiation field can be 
approximated by the expression 

Lad{e)TT 



where Rblr is the effective radius of the BLR, and tt is the Thomson depth of the BLR.^ The 
energy density of the external radiation field in the jet frame is u'^^^ = T'^Uext, where Uext = 
J^deuexti^)- 

The quantities Lad, tt, and Rblr which define the energy density of isotropic external field 
component in equation (9) are not well-known for specific sources. The quantity u'^^^ can also be 



^Note that the factor 2n in the denominator of equation (9), instead of the commonly used 4n, corresponds to 
the fact that we are considering the photon density not outside, but inside the emission region (Atoyan & Aharonian 
1996). 
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estimated from the ratio of the uF^, peak flux of the 7-ray Compton component to the peak 
flux fi^ observed in the lower frequency nonthermal synchrotron component. Assuming that the 
observed 100 MeV - GcV fluxes as observed with EGRET from FSRQs in the flaring state are 
primarily due to Compton-scattered radiation from external photon fields, then 



K.t = 9u'B^, (10) 

Js 



where u'^ = B'^/Stt is the comoving magnetic-field energy density, and 5 > 1 is a parameter that 
corrects for the Klein-Nishina effect in the observed 100 MeV - GeV fluxes.^ For definiteness, 
note that the Klein-Nishina transition occurs when the parameter b = 47eg /(mgC^) exceeds unity 
(Blumcnthal and Gould 1970). This implies that 20 cV UV target photons will be scattered in 
the Klein-Nishina regime in order to produce gamma rays with energies exceeding w 800/(1 -|- z) 
MeV (independent of F). As shown by our numerical calculations, this requires that 5 ~ 3 for the 
parameters of 3C 279 that we consider. 

The magnetic field B' in the blob can be determined by introducing the equipartition parameter 
rj = u'^{l + kpe)/u'^ for the ratio of relativistic electron to magnetic-field energy densities in the jet, 
with the factor kpg = u'^/u'^ correcting for the contribution of nonthermal hadrons to the particle 
energy density. We assume kp^ = 1 in the estimate of B' . The energy density in nonthermal electrons 
is inferred from the measured synchrotron flux density Fg{y) oc z/~" in the range I'o < < 1^1, with 
a = 0.5. This gives the magnetic field 

B-(Gau.s) ^ 130 '^^'^-"1" ^ ''^('^'M^J + ->'^' . („) 

and the equipartition magnetic field corresponds to the parameter rj = 1. Here U13 = i//10^^ Hz is 
the frequency where the flux /s(e) = 10~^°/_io ergs cm~^ is measured. 

An alternative estimate of B' can be derived from the ratio of the synchrotron self-Compton 
(SSC) 7-ray and synchrotron peak fluxes fssc/fs^ — '^'s/'^'b (Sikora 1997; Tavecchio et al. 1998). 
Thus u'b = Au', = Af!\2dl/R'^cS^), where A = {ff/ffsc)- This gives 



This equation assumes that Klein-Nishina effects for the SSC 7 rays are unimportant, because 
the Compton process occurs on softer (radio-optical) photons as compared to the UV photons for 
external Compton scattering. In both equations (11) and (12), the strongest dependence of B' is 
through the unknown Doppler factor. By equating these two equations, we can derive an estimate 



^This estimate does not take into account the different beaming factors for the two processes (Dermer, Sturner, 
& Schlickeiser 1997), but is valid when V ^5. 
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for 6, given by 

^3/8 (. , Nl/4 ^3/16 .7/16 „l/4 ,1/8 

• [(l + A;pe)ln(i^iM)]VMW(d)]V8 • ^ ^ 



2.3. Energy Distribution of Protons and of Escaping Neutrons 

Relativistic neutrons, being essentially unaffected by magnetic fields, can escape from a blob 
if their Lorentz factor 7„ is sufficiently large so that they do not decay into protons on the decay 
timescale T^iju) = To In, where tq ^ 910 s is the decay time of a neutron at rest. Moreover, the 
neutrons must not be converted back to protons in further photomeson interactions with ambient 
photons inside the blob during the crossing time fcross — R'/c needed, on average, for a neutron 
to leave the blob. This leakage must be taken into account in the kinetic equation describing the 
evolution of the proton energy distribution Np{E, t) in the blob at time t (allowing for the possibility 
of a non-stationary source). We develop a formalism for the evolution of the proton distribution 
that properly treats the possible escape of protons while in the neutron state. 

The distribution Np{E,t) can be calculated in the framework of the approximation of con- 
tinuous energy losses. Even though inelasticities up to Kpj = 0.5 can be reached in photomeson 
interactions, the continuous energy-loss approximation is reasonable as follows from the analytic 
solution for the kinetic equation of relativistic particles when both continuous and catastrophic 
energy losses are taken into account (Atoyan 1992b). The kinetic equation for Np{E,t) is then 
given by 

where Qp = Qp{E,t) is the particle injection rate, and Pp^{E) = —{dE /dt)py is the proton energy 
loss rate due to p7 interactions. Note that Pp^ includes both the channels with and without conver- 
sion of the incident proton to the neutron, namely p+j ^ p + ir^ + X and p+j —>■ n + Tr~^ + X. This 
can be done because in photomeson interactions the energy loss rate of a neutron that is produced 
in the latter interaction chain is practically the same as of the proton in the first chain. Therefore 
once the nucleon returns to the proton pool inside the blob, it does not 'remember' that during 
some time At < tcross it was in the 'neutron' state. In other words, the energy loss of protons 
remaining in the blob would depend only on the overall time they spend in the blob after their 
injection. 

The parameter Tcsc = TcsdE) corresponds to the characteristic escape time of protons from 
the blob. This should generally take into account the escape of protons both due to direct escape 
from the blob while in the neutron state, as well as due to diffusive escape. In this section, we 
neglect diffuse escape effects, although this process is taken into account to determine the spectra 
of protons at the highest energies in the numerical calculations. The solution to equation (14) is 

Np{E,t) = -\-f P{et,)Q{et„ti)e^p(- f ^^)dti. (15) 



- 11 - 



The parameter et^ = e{E,t,ti) corresponds to the initial energy of a particle at an earlier time ti 
which cools down to a given energy E by the time t; thus £{E, t, ti = t) = E. It is determined from 
the equation 

r*i dE ^^g^ 



Note that the running time ti can be formally both smaller and larger than the time t at which 
the particle energy E is given, i.e., e is simply a trajectory of such particle in the energy space. 

The characteristic escape time of protons can be evaluated if we recall that ATVes^ = (Np/resc) x 
At is the number of particles escaping the soTirce during the time interval At. It can be written 
also as AiVcsc/At = (,cscNn = ^csc^pni^p-yiE) Np{E,t), where = C (- 0.5) is the probability of 
conversion of a proton into the neutron per one p'j collision, is the frequency of such collisions, 
and ^esc is the probability that the produced neutron would escape the blob prior to its conversion 
back to the proton state. Thus 

^esciE)^pnl^jniE) . (17) 



Tesc(^) 

In order to calculate ^esci consider the solution to the equation 

-Qf = -Q^ i^nnPn-r ^n) - {^npt^n-r + V^d) + Qn (18) 

for the energy distribution of neutrons Nn{E,t) in the blob, assuming instantaneous injection at 
time to = 0) which corresponds to the source function Q^[E, t) = NQ{E)5{t — to). In this equation, 
P„^ and Un"/ represent the overall photomeson energy loss and collision rates of the neutrons, 
respectively, which are approximately the same as the respective rates for the proton photomeson 
interactions. Unlike in eq. (14) for the protons, however, we now have to take into account that 
after each photomeson collision, a neutron either remains neutron with probability (= 0.5) 
while losing energy, or is converted to a proton and returns to the proton pool with probability 
^np = l~Cnn- The ueutron escape term also takes into account the decay of neutrons n — > p+e^ + v 
on timescale = Td{E), with En = Ep. 

The solution to equation (18) is given in general by equation (15). For the instantaneous 
injection function Qn{E,t) = No{E)6{t), equation (15) reduces to 



Nn{Ei,ti) = — ' , . exp 

Ppi\Ei) 



"J" ^np^wy 1 dx 



(19) 



Here we use E\ for the running energy at ti , while reserving E for the initial neutron energy. Note 
that in the continuous energy loss-approximation, E should be the same as the energy of the parent 
proton. 

Because the neutrons escaping the source are the ones which survive during the time ti > 
Across — R'/c aftev injection, the escape probability ^esc is formally the ratio Nn{Ei,ti)AEi/NQ{E)AE 
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at time ti = tcross 
tervals AEi/AE 



Taking also into account that, using equation (16), the ratio of the energy in- 
P{Ei)/P{E), we find that 



exp 



icros 



(20) 



This defines the effective time of proton escape Tesc{E) via equation (17). The evolution of neutrons 
outside the blob is given by equation (18) with the source function Qn = Qn{E, i) equal to the rate 
of escaping protons from the blob in equation (14), i.e., Qn{E,t) = Np{Ep,t)/[{l — Kpj)Tcsc{Ep)], 
where Ep = E/{1 — Kp^) is the mean energy of the protons producing neutrons with energy E. 



3. Results of Calculations 

We calculate the spectra of protons, neutrons, and secondaries from photomcson interactions, 
taking inner jet parameters derived from quasi-simultaneous observations (Wehrle ct al. 1998) of 
the 1996 flare from 3C 279 as typical for FSRQs. For calculations of the photomeson production 
effects in low luminosity blazars, we use parameters of the blob derived from flaring states of the 
BL Lac object Mrk 501 (Catanese et al. 1997; Plan et al. 1998). 



3.1. Photopion Production in FSRQs 

The redshift of 3C 279 is z = 0.538, which implies a luminosity distance c^l — 1-05 x 10^^ cm for 
an = 0-3, = 0.7 cosmology with a Hubble constant of 65 km s^^ Mpc^^. Wc use parameters 
derived for the flaring state of 3C 279 during 1996 February 4-6 (Wehrle et al. 1998). Note that 
parameters appropriate to the extended 3 week quiescent state of the same observing period were 
used earelier in Atoyan &; Dermer (2001). The measured variability time scale t^ar — 1 day for the 
1996 flare of 3C 279 implies a characteristic source size 

R' ^ ct^arS/il + 2) = 1.7 X 10^^W(d)(VlO) cm . (21) 

For calculations of n^(e') = Ug{e')/e' in equation (7), we approximate the flux density Fs{e) = 
/s(e)/e oc observed during the flare of 1996 (Wehrle et al. 1998) in the form of a continuous 
broken power-law function, with indices ai = 0.5 at frequencies fo — 10^^ Hz < v = e/h < 
vi = 10^'^ Hz, Q!2 = 1.45 Sit ui < V < U2 = 10^^ Hz, and 03 = 0.6 at > V2- The uFy synchrotron 
radiation flux reaches its maximum value fs^ = fsi^pk) — 1.7 x lO"^'^ ergcm"^ s~^ at e = Cpk = hvi. 
The flaring flux of the Compton 7-ray component peaks at ~ 500 MeV, and is ~ 15 larger than 
ff'' during the flare (Wehrle et al. 1998). 

Spectral models for 3C 279 show that a complete spectral fit requires that the ratio of syn- 
chrotron to SSC peak fluxes A 2± 1 (Bottcher 1999; Hartman et al. 2001). This implies from 
equation (13) that (5 c± 6, which is just above the lower limit 6 > 5 required from the condition 
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that the gamma rays in the jet be transparent to 77 pair attenuation (Atoyan k, Dermer 2001). 

In the calculations, we treat a range of Doppler factors 5 = 6, 10 and 15, noting that there can 
be uncertainties in our knowledge of A, as well as departures from equipartition. Note, however, 
that the range of the Doppler factors cannot be very large due to the weak dependence of 5 on the 
uncertain parameters. The magnetic field is then derived from equation (11). 

During the 1996 February 4-6 flare, a mean apparent isotropic '47r' 7-ray luminosity of 10^^ 
ergs s~^ was measured with EGRET. Assuming that the power injected in relativistic protons is 
equal to the power observed in 7 rays, we take for the comoving proton luminosity L'^ = 10^^/5^ 
ergs s~^. The injection of protons takes place during the time tinj = 2 days in the observer 
frame, corresponding to the time Stinj/{1 + z) in the comoving blob frame. For the spectrum of 
protons, we take Qp{Ep) oc Ep'^ e'x.p{—Ep/ Emax), where E^ax is defined by the condition that the 
proton gyroradius does not exceed R' and that the proton cooling rate t^{Ep) does not exceed the 
maximum acceleration rate given by eBc/Ep. This acceleration rate is consistent with the maximum 
energy of particles accelerated by relativistic shocks (Achterberg, Gallant, Kirk, &; Guthmann 2001). 
The BLR radius Rblr is deduced from equation (9) with the assumption tt = 0.1. Equation 
(10) is used to estimate the external radiation density 7ig^j(e') in the blob frame, noting that 
€'u'^^f{e') = T'^eUexti^) (e.g., Dermer & Schlickeiscr 2002). This approach results in a dependence 
of the BLR radius on 6 (which we assume to be equal to F), Rblr = 0.11 pc and 0.47 pc for 
6 = 6 and S = 10, respectively. For the range of Doppler factors used in our calculations, the blob 
remains within the BLR during the entire time of injection of relativistic protons. 

Figure 1 shows the spectra of protons and escaping neutrons and gamma rays for (5 = 6 and 
(5 = 10 in Figures la and lb, respectively. The dotted curves show the overall spectra of protons 
injected into the blob, and the thin solid curves show the spectra of protons which remain in the 
blob by the time that the blob reaches the edge of the BLR at stationary- frame time t = Rblr/c- 
The thick solid curves show the overall spectra of neutrons which escape from the relativistic blob 
during the time of its propagation through the BLR. The feature in the neutron spectrum between 
10^^ and 10^^ eV in Figure la and, to a lesser extent, in Figure lb, is due to the neutron energy 
losses on the external radiation field which are only significant above ~ 10^^ eV. The magnetic 
fields for the 5 = 6 and 6 = 10 cases in Figures la and lb are equal to 18.1 G and 7G, respectively. 
The corresponding characteristic maximum energies of the accelerated protons in these two cases 
are equal to 2.5 x 10^^ eV and 6.5 x 10^^ eV in the blob frame, which are translated to energies by 
a factor of F ~ 5 times higher in the stationary frame. 

Neutrons escaping from the blob moving with Lorentz factor F will form a highly coUimated 
beam with opening angle ^ ~ 1/F. These neutrons continue to interact with the external radiation 
field through photopion production, and are also lost through decay. Both effects modify the neu- 
tron spectra until the neutrons reach the edge of the BLR with the external UV field. The spectra 
of surviving neutrons reaching this edge are shown by the dashed curves. Photomeson interactions 
of neutrons outside the blob result in the production of a beam of ultra-high energy neutrinos, 
electrons, and 7-rays with the same opening angle ^ 1/F. The expected number of neutrinos to be 
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detected with ground-based neutrino telescopes is calculated below. The electrons may be deflected 

by the magnetic field in the BLR to produce synchrotron or Compton-scattered radiation, but we 
do not consider this emission further here because the fluxes will depend significantly on unknown 
properties of the medium outside the jet, such as the magnetic field strength and geometry. 

The spectra of vr'^-decay gamma rays produced by neutrons outside the blob which escape from 
the BLR without undergoing subsequent 77 interactions are shown by the dot-dashed curves in 
Figure 1. These gamma rays have very high energies, ranging from lO-*^^ to lO-*^^ eV. The thin curves 
in Figure 2 show the 77 pair-production opacity due to interactions of gamma rays with photons 
of the BLR, calculated for a gamma ray traversing the radius of the BLR. Also shown by the heavy 
curves are the 77 opacities within the blob. This figure demonstrates that gamma rays produced 
within the blob are subject to very strong attenuation, and therefore have no chance to escape. 
The attenuation is also very strong for photons produced outside the blob in the BLR unless these 
photons have very high energies (> 10^^ eV) or are produced close to the edge of the BLR. The 
peak feature in the escaping 7-ray spectrum between 10^^ and 10^^ eV in Figure la is due to the 
threshold character of photomeson interactions of neutrons with the external radiation field, which 
continue to produce gamma rays until they reach the edge of the BLR. 

We calculate the fraction of energy initially injected in the form of protons above 10^^ eV, 
Wp{> lO^^eV), which is taken away from the blob in the form of neutrons. For the cases of (5 = 6 
and 10, as shown in Figure 1, these fractions are 14.8% and 25.9%, respectively. Neutrons with 
energies greater than 10^'' eV that escape from the the BLR will decay on length scales 

Id « 1 (EJIO^^) kpc (22) 

The fraction of injected relativistic proton energy Wp{> lO^^eV) carried by these neutrons, which 
we denote by Wn{> 10^'' eV), is equal to 0.00059% and 3% for 6 = 6 and 10, respectively, for the 
parameters used in Figure 1. The fraction of ultrarelativistic energy that escapes from the BLR in 
the form of gamma rays above 10^^ eV is equal to 0.23% and 5.2% for 6 = 6 and 10, respectively. 
The fact that a smaller fraction of initial energy is taken away by neutrons for (5 = 6 than for S = 10 
is explained by the stronger external radiation field in the former case. This is a consequence of the 
strong dependence of the inferred external radiation field energy density Ug^.^ oc (5~^^/^, as follows 
from equations (10) and (11), which strongly attenuates the beam of neutrons escaping the blob 
close to the black hole. Our calculations show that the fraction of injected energy in ultrarelativistic 
protons, Wp{> lO^^eV), that escape the BLR in the form of neutrons and gamma rays can reach 
~ 10% for different model parameters, when attenuation of the neutrons escaping the BLR is not 
very strong. In particular, the yield in neutrons leaving the BLR is always high if the proton 
acceleration in the blob occurs not during a short flaring period, as in Figure 1, but during the 
time t tsLR = Rblr/c when the blob is traveling through the isotropic external radiation field. 
On the other hand, a higher rate of degradation of neutrons results in higher fiuxes of neutrinos 
produced both inside and outside the blob. Figure 3 shows the fluences of neutrinos integrated 
over several days in the observer frame determined by the time for the blob to pass through the 
BLR, for the same parameters as used in Figure 1. The solid and dashed curves show the fluences 
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calculated for 5 = 6 and 10, respectively. The thick and thin curves represent the fluences of 
neutrinos produced by photopion interactions inside and outside the blob, respectively For the 
spectral fluences shown in Figure 3, the total number Ny of neutrinos that could be detected by a 
1 km^ detector such as IceCube, calculated using neutrino detection efficiencies given by Gaisser, 
Halzen, & Stanev (1995), are 0.29 and 0.078, for (5 = 6, and 0.13 and 0.076, for 5 = 10, where the 
pair of numbers refer to neutrinos formed inside and outside the blob, respectively. 

The dot-dashed and the triple-dot-dashed curves represent the neutrino fluences that would 
be expected for the same time scales and source parameters, but neglecting interactions with the 
external radiation field. The numbers of neutrinos that are predicted to be detected by a 1 km^ 
detector in this case are more than an order-of-magnitude lower, and are 0.030 and 0.0047 for 
6 = 6 and 10, respectively. As previously noted (Atoyan &: Dermer 2001), the presence of external 
radiation fields is crucial for prospects to detect neutrinos from blazars with high-energy neutrino 
telescopes. 

In order to understand how the production of neutrons will change when the inner jet expands, 
we show in Figures 4a and 4b the spectra of protons and escaping neutrons and gamma rays, as in 
Figure 1, but calculated for different sizes of the blob corresponding to tyar = 0.3 and 3 days, with 
S = 10. Here we have fixed the external radiation field density and size of the BLR {Rbrl = 0.41 
pc) to the values derived in Figure lb for 6 = 10. Moreover, we assume that the acceleration 
rate of protons is such that the same amount of total energy as in Figure 1 is injected in the 
protons continuously during the entire period that the jet travels through the BLR. The fractional 
energy taken away by neutrons from the compact blob and from the BLR is 30% and 7.5%, 
respectively, for t^ar = 0-3 d. In the case of the infiated blob with t^ar = 3 d, these fractions drop 
down correspondingly to 12% and 2.4%. 

In Figure 5 we show the radiation fluxes produced in and escaping from the blob following the 
electromagnetic cascade initiated by energetic pion-decay electrons and gamma rays for the case S = 
10. Figure 5a shows the cascade radiation spectra at the maximum, which is reached at time t = 2d 
when the injection of newly accelerated protons in the blob stops. Figure 5b shows the radiation 
fluxes at the time when the blob is leaving the external radiation field region. The thick and thin 
curves correspond to synchrotron and Compton-scattered radiation, respectively. The radiation of 
the first generation of electrons, which includes both the electrons from tt^ decays and the electrons 
produced by absorption of 7r°-decay gamma rays in the blob, are shown by the solid curves. The 
synchrotron radiation strongly dominates the Compton-scattered emission in the first generation 
because the electrons of the first generation are ultra-relativistic, such that Compton scattering 
is in the extreme Klein-Nishina regime. The second, third, and fourth generations of the cascade 
radiation are shown by the dashed, dot-dashed, and triple-dot dashed curves, respectively. In the 
calculations of the escaping cascade radiation spectrum, we have not considered the emission from 
electrons which are produced outside the blob. These electrons become isotropized and produce 
much weaker observed fluxes than the Doppler-boosted emission produced in the blob. Most of 
the r-y-y opacity is due to the external radiation fleld above several GeV (see Figure 2). Thus the 
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spectra of different generations of cascade radiation escaping the BLR are cut off at nearly the 
same energy as long as the blob is in the BLR. At the same time, the cascade radiation made when 
the blob is outside the BLR can extend to higher energies. 

The heavy dots show the total spectrum. 

As can be seen in Figure 5, emission near 1 GeV can escape from the blob and BLR. These 

figures show that in both cases the main contribution to the resulting gamma-ray flux is due to 
synchrotron radiation of the first generation of the cascade electrons born mostly with energies 
above lO^"^ eV. 

Noticeable gamma-ray fluxes may also be produced by secondaries from p7 interactions in the 
blob at later stages when the blob is still within the BLR, as seen in Figure 5b. One may also 
expect that ~ 10 GeV gamma-ray photons produced at these later stages can be due to the escaping 
neutrons interacting with the UV photons outside of the BLR. 

3.2. Photopion Production in BL Lac Objects 

BL Lac objects have generally lower luminosity than FSRQs and have weak or absent line 
emission, suggesting that the column density of the gas in the central parsec region is smaller. This 

may be due to a reduction in the fueling rate as the surrounding gas accretes onto the central engine 
(Bottcher &: Dermer 2002; Cavalicrc & D'Elia 2002). If this scenario is correct, then the mean black 
hole masses in BL Lac objects should be even greater than in FSRQs. Consequently, the maximum 
temperature of the disk in a BL Lac object should be significantly lower than in FSRQs (recall 
that Emax oc M^f^^^ {Mc^ / L^dd)^^^ j where L^dd is the Eddington luminosity and M is the mass 
accretion rate; see Dermer Sz Schlickeiser 2002). In our calculations, we take emax = 3eVfor BL 
Lac objects. 

We consider Mrk 501, a high-frequency peaked BL Lac object that has been observed at TeV 

energies, as representative of this class of sources. Mrk 501 has a dshift z = 0.034, so that cIl = 160 
Mpc. Variability has been detected on time scales as short as 1 hour, although variability on time 
scales of order 1 day are also common. We therefore consider both t^ar = 0.1 and 1 day (Quinn et 
al. 1999; Sambruna et al. 2000). For the compact jet with t^ar = 0.1 d, we approximate the flux 
density using the April 16th, 1997 flare data (Plan et al. 1998), with indices ai = 0.5 at frequencies 
uo = 10^^ < u = e/h < ui = 10^^ Hz, and a2 = I at ui < u < 1^2 = 10^° Hz, and assuming a 
very steep cutoff above U2. The maximum uF^, flux is at the level of ~ 6 x 10"^'' ergs cm"^ s~^. For 
the 'quiescent' state with t^ai- = 1 day, which may correspond to a less powerful but continuous jet, 
we approximate the spectral fluxes with the same power-law indices a, but with break frequencies 
vi = 3 X 10^^ Hz and f 2 = 3 x 10^^ Hz, and a maximum vF^ flux at the level of 3 x 10^^^ ergs 

—2 —1 

cm s . 

Figure 6 shows the total energy spectrum of injected protons and protons remaining in the 
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blob at the time when the jet is at a distance R = 0.3 pc, which we assume to be the effective radius 
of the region where there is an isotropic scattered radiation field. In this calculation, we treat the 
tvar = 1 d case, and let 5 = 10. This results in an cquipartition magnetic field Bgq = 0.37 G, and a 
maximum energy of accelerated protons equal to 9 x 10^^ eV in the comoving frame, which is larger 
by a factor T ~ 10 in the stationary frame. For the external field, we assume that the maximum flux 
of the disk radiation cannot exceed the flux of 3 x 10~^^ ergs cm~^ observed in the optical/UV 
range. The intensity of the isotropic radiation field is calculated assuming tt = 0.01. For these 
parameters of external radiation-field, there is essentially no photomeson production by neutrons 
outside the blob, and consequently very low gamma-ray production outside the blob (dot-dashed 
curve) . 

The production of high-energy neutrons inside the blob through photomeson interactions with 
the internal synchrotron photons is still very effective for the highest energy protons with Ep > 
10^^ eV. The fraction Wn{> W^'^eY) /Wp{> lO^^eV) = 0.12, and practically all of this neutron 
energy is carried from the supposed isotropic external radoiation to kpc scales. The maximum 
energy of protons is, however, smaller compared with the S = 10 case for 3C 279 (Figure 1), which 
implies a shorter propagation distance for the ultra-high energy (UHE) neutron beam. The power 
deposited in the neutron beam during the flare is > 4 orders of magnitude smaller than the power 
of neutrons escaping the BLR in Figure 1. Moreover, for larger Doppler factors, 6 ~ 15-45, which 
are implied by the modeling of TeV and X-ray flares (Plan et al. 1998; Krawczynski et al. 2002), the 
power in neutrons would decrease further. Our calculations show that Wn{> lO^^eV) ^ 2.6 x 10^^ 
and 3.6 x 10^^ ergs for S = 10 and 25, respectively. Note that the fraction of energy carried 
by > 10^^ eV gamma-rays produced by neutrons escaping the blob in the isotropic external 
radiation region in Figure 1 is absolutely negligible, W^{> lO^^eV) /Wp{> lO^^eV) « 4.7 x 10"^ 

In Figure 7, we show the fluxes calculated for the case t^ar = 0.1 d, assuming S = 25. Figure 
7a is calculated for B' = Bgq = 0.5 G, corresponding to cquipartition as defined in equation (11), 
whereas Figure 7b assumes B' = 0.2 B^q = 0.1 G. The latter case may be of interest because the 
spectral fits for the 16 April 1997 flare in Mrk 501 seem to require fields well below cquipartition. 
The maximum energies of accelerated protons in the stationary frame are then Ejnax — 7.5 x 10^^ 
eV and E^ax = 1-5 x lO'''^ GeV in Figures 7a and 7b, respectively, so that the effective maximum 
energy of the neutrons escaping the inner jet (and therefore also the weak external radiation field) 
does not significantly exceed 10^^ eV, as seen in Figure 7. This suggests that in BL Lac objects, the 
transport of energy by ultra relativistic neutrons is still possible, but the maximum distances for 
energy transport by the neutrons will be significantly shorter than in FSRQs. Note also that the 
number of neutrinos to be expected by a km-scale neutrino detector in these two cases is extremely 
low, being less than 10""^ in both Figures 7a and 7b. 
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4. Discussion and Summary 

Observations of powerful gamma-ray flares from blazars with EGRET and ground-based gamma- 
ray telescopes have confirmed that radio-loud AGNs accelerate particles to GeV and TeV energies, 
respectively. Particle acceleration to higher energies may occur, but this cannot be established 
through gamma-ray observations because of attenuation by ambient and diffuse extragalactic op- 
tical and infrared radiation fields. Although Icptonic models have been successful in fitting the 
spectra of blazars, models of first-order Fermi acceleration at both nonrelativistic and relativistic 
shocks imply that acceleration proceeds more effectively for protons than for electrons, because 
the Larmor radius of a particle involved into the acceleration process must be larger than the 
shock width, which is more easily satisfied by hadrons (e.g., Gallant 2002). The large proton-to- 
electron ratio in the galactic cosmic rays supports this contention. Thus our underlying assumption 
that protons are injected with at least the same power as inferred from the measured gamma-ray 
luminosity of a blazar seems reasonable. 

Evidence for the acceleration of relativistic protons and heavier nuclei in blazar jets can be 
provided by three lines of evidence as described in the following subsections. First, hadronic ac- 
celeration could be established indirectly by detection of spectral features in the electromagnetic 
radiation produced by secondaries of the inelastic interactions of high-energy hadrons in blazar 
jets. The second and most compelling line of evidence will be provided through the direct detec- 
tion of high-energy neutrinos by km-scale detectors such as IceCube.^ The third line of evidence 
involves consequences of the production of UHE neutrons and gamma-rays which transport signif- 
icant amounts of energy from the central region of an AGN to large distances (Atoyan & Dormer 
2001). This last feature is sufficiently interesting that we devote a separate subsection to applying 
this model to explain observed features of radio-loud AGNs. Finally we summarize, considering 
also the possibility that ultra-high energy cosmic rays (UHECRs) are accelerated in the inner jets 
of blazars. 

4.1. Electromagnetic Radiation from Nonthermal Hadrons 

Discriminating between hadronic and leptonic origins of the gamma radiation from blazars or 
other classes of sources is possible, although not an easy or a straightforward task. If the emission 
has a hadronic origin, then the emerging gamma-ray spectrum results from a cascade induced by 
high-energy secondary leptons and gamma-rays. Photon spectral indices from cascade radiation 
tend to be rather hard, and take values between 1.5 and 2.0 in the medium-energy gamma-ray 
regime explored by EGRET and GLAST ^ (see Figure 5). But this is also a common feature 
for gamma-ray spectra from an electron-photon cascade independently of whether the cascade is 



^http: / /icecube. wisc.edu/ 
*http: / /glast.gsfc. nasa.gov 
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initiated by electrons and gamma-rays of secondary (pion-decay) origin or by primary (directly- 
accelerated) electrons. 

A principal difference between the hadronic and leptonic induced cascades does, however, exist. 
It consists in a very significant contribution of the synchrotron radiation by the ultra-relativistic 
electrons (including the ones from the absorption of 7r*^-dccay gamma-rays) of the first cascade 
generation in the high energy gamma-ray flux in the case of hadronic initiation of the cascade. The 
Lorentz factor of these electrons 7 S> 10^ is much higher than the maximum possible Lorentz factor 
for primary electrons, which is limited by synchrotron losses in the strong magnetic fields S > 0.1- 
10 G characteristic for the inner jets of blazars. As a result, the synchrotron radiation of these 
electrons can greatly exceed the characteristic maximum energy max ~ 25-100 MeV (Guilbert, 
Fabian, & Rees 1983; e.g., de Jager et al. 1996) for synchrotron emission of directly accelerated 
electrons, as is apparent from Figure 5. Additionally taking into account the Doppler boosting of 
the energy £^s,max when transforming to the observer's frame, we can say that confirmation of a 
significant contribution of hard-spectrum synchrotron flux in the gamma-ray flares at > 0.1 GeV 
with GLAST could become a strong argument in favor of UHE hadron acceleration in the jets. 

A synchrotron contribution could be confirmed experimentally by the detection of a significant 
polarization in the total flux, but neither GLAST nor other high energy gamma-ray detectors yet 
have the capability for the relevant measurements. It is possible that the synchrotron origin of the 
GeV radiation could be revealed by accurate modeling of the spectral and temporal behaviors of 
the broad-band gamma ray fluxes detected by GLAST, as well as by forthcoming ground-based 
Cherenkov detectors. Of particular interest in this regard could be the interpretation of flares 
which would decline rapidly at 0.1-1 GeV energies at the flare fading stage. The key observation 
here is that in FSRQs, the Compton radiation at these energies should be due to upscattering of 
external UV photons by electrons with £" < 1 GeV, which could have cooling times larger than 
the measured decline time, whereas the GeV synchrotron radiation is due to electrons with many 
orders of magnitude higher energies, and so it would have much shorter cooling times. 

If radio galaxies are misdirected blazars and steep spectrum radio sources are off-axis FSRQs 
(Begelman, Blandford, &: Rees 1984; Urry & Padovani 1995; Orr &; Browne 1982), then the gamma- 
ray spectra of radio galaxies could also reveal hadronic acceleration in blazar jets. Nonthermal 
cosmic rays that escape from the inner jets of radio galaxies can accumulate in the central gas-rich 
region of the AGN. Nuclear collisions of these cosmic rays with ambient gas would produce a quasi- 
stationary (on time scales of thousands of years) and compact gamma-ray halo with a characteristic 
TT^-decay feature near 70/(l+z) MeV. Note however that the bremsstrahlung emission from tt^- 
decay electrons of ~0.1-1 GeV energies is likely to conceal the 70 MeV benchmark (Schlickeiser 
1982). Nuclear interactions will give secondary 7r-decay emission that is brightest in the central 
BLR core of a radio galaxy where the target particle density is highest, and dimmer in the extended 
region due to the tenuousness of the interstellar medium on scales of hundreds of pc. Because of 
the rather small size of this hadronic halo, the gamma-ray detectors would detect it as a point 
source. Such a compact hadronic halo could be still resolved in the radio band, where the radio 
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fiux produced by secondary electrons of GeV energies would show an increase in brightness towards 

the center of the AGN on the milliarcsecond scale. Note that at very high energies, a halo on the 
much larger scale (multi-Mcgaparsec) could be formed around such sources due to absorption and 
cascading of multi-TeV gamma-rays in the intcrgalactic medium (Aharonian et al. 1994). 

Considering possible signatures of rclativistic hadrons during the flares, we note that a pro- 
nounced 70 MeV Tr'^-decay maximum should be due first of all to hadrons with kinetic energies ~ 1- 
10 GeV; therefore it would imply a mass-loaded jet. As we show in the Appendix A, the mass-loaded 
jet models are rather unlikely because they confront severe energetics difficulties. Nevertheless, the 
hadronic 7r°-decay gamma-ray feature, although significantly broadened and Doppler-shifted, could 
be expected in those particular cases when the inner jet dumps its energy in collisions with the 
dense target, such as a broad line cloud crossing the jet. Observations of radio galaxies with imaging 
X-ray telescopes, GLAST, and ground-based air Cherenkov telescopes with low energy thresholds 
will provide spectral data that should be analyzed for evidence of nonthermal hadrons in light of 
these considerations. 



4.2. Neutrinos from FSRQs and BL Lacs 

Detection of neutrinos from blazars will provide the most compelling evidence for hadronic 
acceleration in blazar jets. For photohadronic jet models, we have argued here and elsewhere 
(Atoyan & Dermer 2001) that the presence of a strong UV accretion-disk radiation field is required 
for blazar neutrino fiuxes to be detectable with km-scale high-energy neutrino telescopes. The disk 
radiation field is needed to produce a quasi-isotropic scattered radiation component produced in 
the BLR of FSRQs. The existence of an external radiation field in FSRQs is suggested by luminous 
broad emission lines detected from these objects. For UV fields expected from optically-thick 
accretion disks, photomeson interactions can efficiently extract the energy of relativistic protons 
with (relatively modest) energies > 10"*^^ eV, resulting in the production of neutrinos with E^, > 
3 X 10^^ eV. For flares like the one detected in February 1996 from 3C 279 with EGRET, the fluence 
of neutrinos shown in Figure 3 is at the level of 10~^ ergcm"^, implying ~ 0.3 per flare (i.e., a 
detection probability of ^ 30%). This would suggest a realistic possibility for IceCube to detect a 
few to several neutrinos from several flares over the course of a year. Given the flaring duty cycle 
of FSRQs estimated from EGRET observations (~ 10-20%), this implies that a few high-energy 
neutrinos will be detected from FSRQs such as 3C 279 on a time scale of 1 year. Considerably 
more will be detected if the efficiency to accelerate hadrons exceeds our baseline assumption given 
by equating proton and gamma-ray power. The additional possibility that some neutrinos will also 
be detected during quiescent states strengthens our prediction that FSRQs will be detectable high- 
energy neutrino sources with km-scale neutrino telescopes. We note that on a 1 year observation 
time scale, detection of only two neutrinos with energies above 30 TeV from the same direction 
pointed to a known blazar would represent a very significant detection of an astrophysical source, 
because at these energies the background cosmic-ray induced neutrino flux expected within one 
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square degree (which is typical for the angular resolution of a high-energy neutrino telescope in ice 
or water) is < 0.03 (Gaisser, Halzen, & Stanev 1995). 

For BL Lac objects, where the accretion-disk and scattered radiation components are much 
weaker, and with GeV gamma-ray fluxes that are also typically not as bright as in FSRQs, the 
chances for detection of neutrinos even with a 1 km-scale neutrino detector are negligible. For the 
proton fluxes shown in Figures 6 and 7 calculated for parameters of the April 1997 flare in Mrk 
501, the probability of neutrino detection during the flare is less than 10~^. It is also important to 
note that unlike the case of p7 interactions with the internal synchrotron photons, the intensity of 
P7 interactions with the external radiation does not depend on the blob size. Therefore in FSRQs 
the neutrino production from photomeson interactions with the external radiation field does not 
decline so much in the process of expansion of the blob as it propagates through the quasi-isotropic 
external radiation field at R < -Rblr as it does in BL Lac objects. We note that the blob expansion 
seems unavoidable because of the difficulty in confining the high internal pressures of the blob by 
the pressure of the external medium. 

An important issue concerns the quenching of the jet as it sweeps up mass from the BLR. 
As BLR clouds execute Keplerian rotational motion about the central supermassive black hole, 
some will pass through the path of the jet outflow. This will cause energization and deceleration 
of the radiating plasma, as in the external shock model of gamma-ray bursts (Meszaros & Rees 
1993). Suppose that the AGN jet is beamed into a fraction fh = 0.01/_2 of the full sky. In a 
simple model with a uniform BLR gas density ublr within a spherical BLR with Thomson depth 
tt = ublrO'tRblr, the total mass in the beaming cone is Mj{< Rblr) — 525O/bTT-Rf8-^0' 
where Ris = Rblr/^0^^ cm (see Blandford Sz Levinson (1995) for BLRs with radially dependent 
densities). The time for BLR material to pass through the cone of the jet with opening half-angle 
0j = 2\/7fe is tj = 1.1 X 10^ f^^2 Rii^ / ^9^"^ where the mass of the central supermassive black 
hole is Mg = Mhh/W^ Mq. Assuming that the directional power of rclativistic plasma expelled in 
the AGN jets is comparable to the apparent isotropic gamma-ray luminosity L^, then the flow will 
be significantly decelerated when L^tj = r^Mj (< Rblr)c^, where F is the Lorentz factor of the 
plasma outflow. Thus, unless the mean apparent isotropic luminosity exceeds 



the jet will be quenched by sweeping up BLR gas. This quenching process could make it difficult 
to form extended relativistic jets on the 10-100 kpc scale, as may be required to explain enhanced 
X-ray emission in the knots of PKS 0637-752 (Tavecchio et al. 2000). The transport of inner jet 
energy to the extended jets via neutral beams avoids this difficulty. The neutrons and gamma- 
rays will pass through the BLR material without attenuation unless the Thomson depth of BLR 
is very high, with tt > 10. In these cases the jet and neutral beam outflow would be quenched, 
the nonthermal radiation would be reprocessed into lower energy emissions, and only the neutrino 
signal would penetrate through the dense inner region without reprocessing. If so, the neutrino 
emission from such buried jets in non-blazar sources could reveal these types of AGN. We return 
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to this speculation in Section 4.4. 

4.3. Energy Transport to the Extended Jet by Neutral Beams 

Relativistic neutrons produced by jrf collisions in the compact inner jet are no longer confined 
by the jet's magnetic field, and so can escape from the blob unless the neutron decays or interacts 
with another soft photon that results in conversion to a proton inside the blob. Our calculations 
show that in FSRQs with luminous external radiation fields, up to ~ 30% of the total energy 
injected into the blob in the form of relativistic protons with E > 10^^ eV can leave the blob in the 
form of relativistic neutrons with Lorentz factors 7„ > 10^. 

For the case of FSRQs, neutrons escaping the inner blob continue to propagate in the BLR 
with its quasi-isotropic scattered accretion-disk radiation up to distances R ~ Rblr ~ 0.1-1 pc. 
Subsequent interactions of the neutron beam with the scattered disk photons will result in the 
production of beamed secondaries. Neutrino secondaries will add to the overall neutrino flux. 
Importantly, a beam of gamma-rays with energies above 10^^ eV is also produced. Because these 
gamma-rays are made outside the blob and propagate radially outward in the blob's frame, they do 
not interact with the internal synchrotron photons, and therefore have a chance to escape from the 
central region of the quasar (see Figure 2) . Depending on the extent and intensity of the isotropic 
external radiation field, as much as several per cent of the accelerated proton energy Wp(> 10^^ eV) 
injected in the inner jet at distances R < Rblr could thus be converted to gamma-rays which escape 
from the BLR. 

Figure 8 shows the mean-free-path of gamma rays due to pair-production interactions with 
photons of the diffuse cosmic microwave background radiation for z = (solid curve) and z = 0.538 
(dashed curve) . This calculation also takes into account the diffuse extragalactic infrared radiation 
field. We approximate this field by a spectrum uJi, oc const at the level of total energy density 
uiR = 0.01 eV cm~^ between 0.001mm and 1mm (e.g., Primack et al. 1999), though the precise 
value is not important on the size scales < 1 Mpc under consideration. 

Our calculations show that more energy is typically carried out from the BLR by relativistic 
neutrons than by gamma rays. This is especially true if the degradation of the neutron beam 
escaping the blob but still in the BLR is not very strong, which will always be the case for neutrons 
escaping the blob at distances R ~ Rblr- The energy carried away by neutrons with En > 10^^ eV 
may reach ~ 10% of Wp{> 10^^ eV) for accelerated protons with number index ap ~ 2. Adding 
the energy of the gamma rays to this, we find that ~ 5% can be used as a rough estimate for the 
fraction of total (starting from 1 GeV) energy of accelerated protons in the inner jet that is carried 
by the collimated beam of UHE neutrons and gamma rays to >kpc distances in FSRQs. 

For a flare like the one in February 1996 from 3C 279, the gamma-ray fluence measured with 
EGRET corresponds to an apparent isotropic energy release of 2 x 10^^ erg. Assuming that the 
injection power in nonthermal protons in the inner jet is the same as for the parent electrons, the 
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overall total energy in the neutral beam emerging from the inner jet per a single such flare could 
reach the level ~ 10^^6^q eig in the stationary frame. This energy will be collimated in the cone 
with opening angle ~ 1/5, and it will be deposited at distances up to 0.1-1 Mpc, provided that the 
spectra of ultra-high energy neutrons extend up to 10^^-10^^ eV, as in Figure 4b. 

After neutron-decay or gamma-ray attenuation, the charged particles emerging in the direction 
of the parent neutral beam will interact with the ambient extragalactic medium via local magnetic 
and photon fields. In the case of the neutrons, a small fraction ~ mg/mp of the neutron energy 
will appear in the form of /?-decay electrons with Lorentz-factors 7 > 10^. This energy will be 
immediately available for radiation in the synchrotron and Compton processes with uF,y peaks 
at hard X-ray and multi-TeV energies, respectively, for a 1 /uG field and a dominant local CMB 
radiation field (Dermer 2002). The dominant fraction of the initial neutron energy remains, however, 
in the protons. 

Synchrotron radiation of ultrarelativistic neutron-decay protons could in principle explain the 
extended Chandra X-ray jet emission if the ambient magnetic field is sufficiently large (Aharonian 
2002). Proton synchrotron radiation does not, however, provide a mechanism for deposition of 
energy of the neutral beam into the surrounding medium to drive shocks and accelerate ~ 1-100 
GeV electrons, which is required to explain the radio and optical emission from the extended jets. 

In order for the transport of energy by the neutrons to be efficient from the point of view of 
powering the extended jets in radio galaxies, there should be a mechanism for effective transfer 
of the neutron-decay proton energy to the surrounding plasma. Such a mechanism is expected 
through interactions of the relativistic protons as well as the e+-e~ pairs from the gamma-ray 
beam with the ambient magnetic field in the surrounding medium. In the case of a significant 
transverse magnetic-field component B±, both the protons and leptons will change their initial 
directions on the Larmor timescales ti = E/eB±c, which implies that a significant fraction of their 
momentum and energy will be transferred to the ambient medium. Another, and possibly even 
more effective way for transferring the initial momentum and energy of the charged particles is 
by generating MHD turbulence through beam instabilities excited by the highly beamed neutron- 
decay protons and charged secondaries emerging from the charge-neutral beam. The continuous 
transfer of momentum and energy from the beam to the intergalactic medium could stretch the 
ambient magnetic field to form a channel facilitating beam propagation, or drive the medium into 
relativistic motion in the forward direction to create a powerful relativistic shock. 

In the latter scenario, relativistic shocks can accelerate ambient electrons (as well as protons), 
thereby injecting a power- law electron spectra Qe{E) oc E""' with a > 2 in the downstream 
region. Synchrotron radiation from these electrons could then explain (Dermer and Atoyan 2002) 
the broad-band nonthermal radio, optical, and X-ray emission from hot spots and knots in the 
extended X-ray jets of extragalactic jet sources observed with the Chandra Observatory. The idea 
that the observed X-ray jets might be due to propagation of a powerful beam of gamma rays 
has also been considered by Neronov et al (2002) (see also Dermer & Li 2001), who suggested 
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a "non-acceleration" scenario, where the observed X-rays could be due to synchrotron radiation 
from electrons formed in the electromagnetic cascade initiated by ultra-high energy gamma rays. 
Our calculations show that although some contribution to the observed X-ray fluxes from the 
cascade electrons is not excluded, the spectral features of the observed radiation, in particular. 
X-ray energy spectral indices ax > 1, might be difficult to explain in a pure cascade scenario. 
Moreover, interpretation of the radio and optical fluxes from the same X-ray knots apparently 
makes unavoidable the requirement of in-situ acceleration of electrons in the X-ray knots. Given 
the total energy released by the neutrons and gamma rays per single flare, the energetics of these 
kpc-scale knots could be explained by superposition of jet activity on time scales up to thousands 
of yrs (since Ikpc translates to ~ 3 x 10'^ yr of light travel time). Thus, the bright knots might 
represent the observable signatures of the past history of activity of the central engine of the blazar 
when powerful acceleration of UHE cosmic rays had occurred in the inner jets of the source. 

Considering now the BL Lac objects, we note that in the case of weak or negligible external 
radiation field components, any significant transport of the inner jet energy by ultra-high energy 
gamma-rays in those objects is practically absent. This is because gamma-rays arc not produced 
through n7 interactions of the neutrons outside the compact blob, whereas the ultra-high energy 
gamma-rays produced in p7 collisions inside the blob will be absorbed due to 77 collisions with 
the same internal synchrotron target photons (recall that a^^ ^ CTpj). Nevertheless, up to ~ 10% 
of the energy injected into the inner jet in relativistic protons with E > 1 PeV could still be taken 
out of the blob by neutrons with E > 10^^ eV due to collisions with synchrotron photons while 
the blob remains compact. Because of the latter limitations, and because BL Lac sources are 
generally much less powerful than FSRQs, the total energy that the neutron beam can extract 
from the accelerated protons is much less for BL Lacs than for FSRQs. Moreover, because of the 
smaller sizes and magnetic fields inferred for BL Lac objects, the level of neutron production and 
the maximum energy of accelerated protons is reduced in the former class of sources. Thus the 
neutral beam power is much less in BL Lac objects than in FSRQs. 



4.4. Explanation of Observed Features of Radio-Loud Active Galactic Nuclei 

As a consequence of processes occurring in the inner jets of blazars, we predict that powerful 
collimated neutral beams that extend as far as hundreds of kpc from the nucleus arc ejected in the 
axial jet directions in high luminosity FSRQs. The linearity and the spatial extent of the cjccta is 
a consequence of the decay properties of the neutrals. This process can explain the origin of X-ray 
jets detected by Chandra, which appear straight on scales of up to Mpc, such as the jet in Pictor A 
(Wilson, Young, and Shopbell 2001). The linearity of the radio jets in Cygnus A (Perley, Dreher, & 
Cowan 1984) may also be a consequence of energetic neutral beams powering the lobes of powerful 
FR II galaxies. This interpretation avoids firehose instabilities in models invoking magnetic fields 
to collimate the extended FR II radio jets, or problems associated with jet quenching, as discussed 
in Section 4.2. 
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Our model predicts that lower luminosity blazars, in particular, the lineless BL Lac objects, 
have orders-of-magnitude weaker neutral beams than FSRQs. The maximum energies of the rcl- 
ativistic neutrons are also about an ordcr-of-magnitude less in Mrk 501-like BL Lac objects than 
in 3C 279-like FSRQs. Consequently the X-ray jets driven by the high-energy neutron beams will 
be generally fainter (despite their relative proximity) and shorter in BL Lac objects. This could 
explain the faint kpc-scale X-ray jets detected from objects like 3C 371 (Pesce et al. 2001), which 
is an off-axis BL Lac object (Miller 1975). 

We adopt the point of view, following Urry k, Padovani (1995), that FR I and FR II galaxies 
are the parent populations of BL Lac objects and FSRQs, respectively. The less powerful beams of 
neutrons that occupy smaller spatial scales in BL Lac objects compared to FSRQs conforms to the 
measured relative power requirements in FR I and FR II galaxies (Fanaroff &: Riley 1974). If FR 
I/BL Lacs have smaller Lorentz factors than FR II/FSRQs, as also inferred from observations of 
these sources (Urry k. Padovani 1995), then the weak, broad twin-jet morphologies in FR I galaxies 
can be qualitatively understood. We point out that even though the inner parsec regions of a BL 
Lac are much less dense than the BLRs in FSRQs in view of their relative emission line strengths, 
the jet quenching problem (Section 4.2) can still be severe for BL Lac objects because of their 
weaker jets. 

The transport of energy by neutral particles to the extended jet will be enhanced by an intense 
external radiation field in the inner regions of the jet. The total power into the jet depends foremost 
on the accretion rate, which depends on the amount of gas and dust available to fuel a supermassive 
black hole. Thus our model is consistent with the evolutionary scenario proposed by Bottcher & 
Dermer (2002) and Cavaliere & D'Elia (2002) to explain the spectral sequence of FSRQs, and low- 
and high-frequency peaked BL Lac objects. According to this picture, FSRQs evolve into BL Lac 
objects due to a reduction in the accretion rate of the surrounding dust and gas that fuels the 
central black hole. (Note, however, that the jet power is derived primarily from mass accretion in 
the scenario of Bottcher & Dermer (2002) and from the spin energy of the black hole in the scenario 
of Cavaliere k. D'Elia (2002).) Maraschi k Tavecchio (2002) argue that this picture accounts for the 
correlation between the radio jet and accretion power, where the latter is inferred from emission line 
strengths. High-energy neutrino observations will be crucial to test this scenario, as neutral beams 
of high energy neutrons and gamma-rays will be made in association with high-energy neutrinos. 
Finally we speculate that if ultraluminous infrared galaxies and quasars evolve into radio-loud 
AGNs (Sanders et al. 1988), then some ultraluminous infrared galaxies could harbor buried jet 
sources that will be revealed through their high-energy neutrino emission. 



4.5. Summary 

Production of narrow beams of ultra-high energy neutrons with E > 10^^ eV becomes possible if 
relativistic compact jets of blazars accelerate cosmic rays to these energies. Such neutral relativistic 
particle beams can be especially powerful in FSRQs, where the neutrons are efficiently produced 
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in p7 interactions with photons of the external quasi-isotropic radiation field on timescale of years 
(since the broad-line region size scale Rblr ~ 0.1-1 pc). In radio-loud jet sources, the beam could 
contain a significant fraction of ultra-high energy gamma-rays produced in the BLR through 
photomeson interactions by neutrons which escape from the compact relativistic plasma blob that 
forms the inner jet. These charge-neutral beams of UHE neutrons and gamma-rays (including also 
neutrinos which, however, do not contribute to the process of energy transport) do not interact 
with the extragalactic medium until either the neutrons decay or the gamma rays collide with soft 
photons on scales from kpc to Mpc. The main soft-photon target for 77 collisions is provided by 
the CMB radiation and, in some sufficiently powerful sources, also by low-frequency synchrotron 
radiation self-produced along the jet (Neronov et al 2002)). 

The weaker neutral beams in FR I radio galaxies and BL Lac objects, compared to FR II 
radio galaxies and FSRQs, could explain the morphological differences between the two classes of 
radio galaxies. This scenario assumes that the inner jets of blazars are effective accelerators of 
high-energy cosmic rays. The limits on acceleration that we consider do not preclude the accel- 
eration of cosmic rays to energies > lO^'^ eV in blazars, especially in the more luminous FSRQs. 
This mechanism of relativistic energy transport from the central engine of AGN to large distances 
avoids problems of jet quenching in the central gas-rich regions of AGN, problems of adiabatic and 
radiative losses, and problems of long-term confinement and stability (on scales up to 10^ yr) of 
the jet against different plasma processes. Beam instabilities can, however, provide an effective 
method to transfer the energy of ultra-relativistic protons emerging from the neutron beam into 
the surrounding medium, and driving relativistic shock which ultimately accelerates nonthermal 
electrons to radiate observable synchrotron and Compton emission. Given also the large amount 
of energy that such a beam may contain, which can be of the order of a few percent of the inner 
jet power, our results suggest that the detection by the Chandra Observatory of the narrow X-ray 
jets that remain straight on scales beyond several hundreds of kpc provides observational evidence 
that the relativistic inner jets of blazars, first of all the flat spectrum radio quasars, are powerful 
accelerators of UHE cosmic rays. Detection of high-energy neutrinos from FSRQs with IceCube or 
a Northern Hemisphere high-energy neutrino telescope will provide strong support for this scenario. 

We thank Hui Li for discussions about the physics of jets, and the referee for a very useful 
report. A A appreciates the hospitality and support of the NRL High Energy Space Environment 
Branch during his visit when this work has been done. The work of CD is supported by the Office 
of Naval Research and NASA grant No. DPR S-13756G. 

A. Comparison of Mass-loaded and Photopion Hadronic Jet Models 

Relativistic neutrons are produced with ~ 50% probability in the inelastic collisions of relativis- 
tic protons either with an ambient thermal proton or with a photon above the photopion production 
threshold. However, the mass-loaded jet models, i.e. the ones invoking pp nuclear collisions, can be 
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practically excluded on the basis of efficiency arguments. As equation (1) shows, nuclear interac- 
tions could provide energy loss time scales shorter than the dynamical time t^^^ = R' /c only if the 
source is very thick to Thomson scattering, Tsc > 40. If one requires that, say, non-thermal X-ray 
flux are produced during the flare, i.e. Tsc< 1, then the observed fast decline of the fluxes could be 
attributed to the adiabatic losses of proton energy and fast decline of the gas density on dynamical 
time scale due to rapid expansion of the source. However, the overall energy losses of protons in 
nuclear interactions during that time cannot exceed ~ 2.5%. In practice this means that the real 
power injected into the blob in the form of relativistic protons should be by a factor of > 100 larger 
than the power inferred from the observed radiation fluence. 

Besides these inefficiency problems, a serious argument against hadronic nuclear interaction 
models is that the mass-loaded jets would require too great a total energy in the blob including 
the kinetic energy. Thus, for a blob size R' ^ t^arcS/{l + z) with typical t^ar — Id and (5 ~ 10 the 
Thomson scattering depth Tgc = 1 implies a gas density (in terms of 'H-atoms', and omitting here 
the dependence on the source redshift z) nn ~ 6 x lO^cm"^. The jet mass then is Mj 2± 3.5 M©, 
and the kinetic energy Wkin = (F — l)Mj(? 6 x 10^^ erg for F ~ 5 10. Thus, even such a 
jet with reduced efficiency should have a significant mass load leading to an unacceptably high 
energetics required in a single blob. It would be formally possible to bring the high kinetic energy 
requirements alone (e.g., Aharonian et al. 1999) down to a formally acceptable values by reducing 
further the density of the gas. However, then the total energy of relativistic protons Wp in the blob 
should be increased in order to provide a given level of radiation fluxes by pp-inter actions. 

Indeed, the integral luminosity of the radiation that could be produced in pp-inter actions, 
is given by L'^. ~ VrWp/tpp, where rj is the fraction of relativistic proton energy that goes to 
secondary electrons and gamma-rays. Assuming that most of the electron energy will be processed 

to electromagnetic radiation, we can have rj = 0.5, leaving the rest for the neutrinos. Transforming 
from the blob frame to the observer (or stationary) frame, the total energy Wtot = Wp + W^in can 
be written as 

- 4dVtot(l + ^)\ .,-1 , {U^rCSf , 

where /tot is the total energy flux, i.e. J F{y^diV which is typically larger than the vF{y) that a 
model is intended to explain. Here we have assumed for the kinetic energy (F — 1)/F ~ 1, and 
substituted Kpp = 0.5 for the inelasticity. The energy VFtot reaches a minimum at the density 

~ 2.25 X 10^ ^28/^10 t^^^^^d) 6-'/^ (1 + zf/^ cm'^ , (A2) 

where /_io = /tot/lO"''^^ ergcm"^ s"-*^, and ^28 = cm. The minimum of Wtot is then 

Wrrnn ^ 1-6 X lO^^ d28f'Jio *v4' (d) ^'w + ^)"'^' ^rg . (A3) 

Here 6io = 6/10, and we have assumed F = 5 in equation (Al). Interestingly, VFmm is reached when 
the kinetic and internal CR energies become equal. Note that any noticeable deviation of the real 
density in the jet from nn.min required for Wmin will immediately increase further the total energy 
requirement given by equation (Al). 



Wtot - 47rFmpc2 
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A hadronic nuclear interaction jet model which could make any significant contribution to 
the fluxes of flares observed by EGRET from FSRQs, would have to assume an absolute (in the 
stationary frame) energy in a single flare exceeding by orders of magnitude an energy of a supernova 
explosion. Even for a less powerful but much closer blazars, like BL Lac objects Mrk 421 or Mrk 
501, any significant contribution to the total radiation flux from mass-loaded nuclear interaction jet 
models can be practically excluded. Indeed, even taking into account that the variability timescale 
might be as short as 1 h instead of 1 day, for interpretation of the maximum fluxes at the level of 
fe ~ 10~^ ergcm~^s~^ the minimum energy Wmin would still significantly exceed the 'supernova' 
level 10^"*^ erg, which is problematic to accept for these objects with a relatively modest average 
bolometric luminosities. The only type of relativistic extragalactic sources for which the hadronic 
nuclear interaction models could still be reasonable from the general total-energy consideration 
seem GRB sources where variability timescale is very short (hence the densities could be very high) 
although the peak fluxes arc much higher than for relativistic jets. Noticing that /tot x t^ar is 
actually a good estimate for the fluence $ observed during the flare, we can re-write equation (A3) 
as 

WVin ^ 6.3 X 10^9 das^L^' (tvar/ls) (5/100)'/^ (1 + z)-'/^ erg, 

where $-5 = <I>/10~^ ergcm"^. Obviously, this is a reasonable value for GRB source, although 
the question on the efficiency of the pp-interaction model in case of simultaneous detection of 
non-thermal X-rays and gamma-rays should be still addressed. 



B. Relative Importance of the Direct Disk and Scattered Radiation Fields for 

Photomeson Production 

Equations (8) and (9) can now be used to estimate the relative importance of the photomeson 

interactions of relativistic miclcons with the direct and scattered components of the accretion disk 
radiation, respectively. For the direct disk radiation photon spectral density n(e, Q) = u(e, $7; /i)/e, 
we analytically estimate the dimensionless number of photomeson collisions per unit path dh (in 
terms of gravitational radius Rg), Kdir = i^pjRg/c, where the collision rate Vpj is given by an integral 
(2) without the inelasticity term Kp^. This yields the simple formula 

Kdirih, E) ^ 4.1 X 10^ la/-^^^ (tT^) ' ■ (Bl) 

^"^^ ' ' "■'^ /i,9/4 VlOOeV/ 200 /xb ^ ^ 

Here lad = Lad/LEdd is the total luminosity of the disk radiation in units of Eddington luminosity 
LEdd-, and 

f 1 — cos ih 

Hh,E)= J —^dcosV, (B2) 

where the integration is to be performed only over the collision angles i/j satisfying the threshold 
condition 
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with 7p^„ = E/mpC? (for m„ = mp) corresponding to the Lorenz-factor of the incident proton or 
neutron. For simpHcity of the analytic estimate (Bl), we here approximate the fry interaction cross 
section o"(ei.) as a single step-function with a height a ^ 200 /xb between a\ and (T2, starting from 
the threshold energy €th — 150 MeV. 

A similar analytic estimate for the coefficient Kiso{E) is derived for the scattered quasi- isotropic 
component of the external radiation field (9). For distances h < vblr, i-e. inside the BLR with 
dimensionless radius tblr, this results in 

.,UE) 1.1 X 10 (^) ^ , (B4) 

where emax = ^*/'^^^^ ■> with > 6, as described above. 

Because the scattered radiation is quasi-isotropic, the interaction threshold condition for pho- 
topion production in this field is the same at all heights h < rsLR, and is defined by head-on 
collisions of protons and neutrons with photons with maximiini energy 6 — ^max • For ^.^nax 

~ 20 eV, 

this requires nucleon energies E > Eth — 4 x 10^ nipC^. Meanwhile, because of the follow-on col- 
lisions that nucleons, which are moving predominantly along the jet away from the black hole, 
make with the photons of the direct disk radiation field, the interaction threshold of the nucle- 
ons with this field significantly increases with increasing height h. From equation (B3), we find 
E > Efiir{h) ~ 4.5 X lO^/i.^/'^rripC^, taking into account that the function reaches a maximum 
= 0.331 at ^ = 60°. This means that only particles with energies E ^ 10^^ eV have effective 
interactions with the direct disk radiation component once the jet reaches heights h > 100. Mean- 
while, interactions of all particles with energies above 4 PeV can effectively continue until heights 
h ~ rBLR ~ 2 X 10^-10^ (corresponding to Rblr ~ 0.1-1 pc for a black hole mass Mbh ~ 10^ Mq). 

The characteristic photomeson opacity Tiso — nisofBLR due to interactions with BLR photons 
can be very significant for l^a ^0.1 and tt ~ 0.1. For comparison, the characteristic photomeson 
opacity T^ir = i^dir h from the direct disk radiation field on scales h can be significant (> 1) only 
rather close to the black hole, at h < 100. This applies only to particles well above the interaction 
threshold Edirih), when the integration limits in the integral (B2) cover a significant fraction of 
the maximum allowed range < cos^ < 1, when I{h,E) reaches a maximum about 0.25. Given 
these considerations, we consider only interactions with the quasi-isotropic (i.e, scattered BLR field) 
component of the external disk radiation. Note that the acceleration of nonthermal particles may 
occur at large distances from the central black hole if due either to collision of relativistic shells or 
to interactions of the outflowing plasma with ambient gas in the BLR. 
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Fig. 1. — Spectra of protons injected into the blob (dotted curves), protons which remain in the 
blob when it reaches the edge of the BLR (thin solid curves), neutrons escaping from the blob 
(thick solid curves), and escaping neutrons (dashed curves) and gamma rays (dot-dashed curves) 
which reach the edge of the BLR for the cases of 5 = 6 and (5 = 10 in the upper panel (a) and lower 
panel (b), respectively. Model parameters relevant to the February 4-6 flare from 3C 279 detected 
by EGRET are used, as described in Section 3.1. 
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Fig. 2. — Pair production opacity due to 77 interactions inside the blob and in the BLR, as 
indicated in the figure legend. 
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Fig. 3. — Fluences of neutrinos integrated over several days in the observer frame determined by 
the time for the blob to pass through the BLR, for the same parameters as used in Figure 1. The 
solid and dashed curves show the fluences calculated for (5 = 6 and 10, and the thick and thin curves 
represent the fluences of neutrinos produced by photopion interactions inside and outside the blob, 
respectively. The dot-dashed and 3-dot - dashed curves show the fluences due to jrf collisions if 
external radiation field is not taken into account. 
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Fig. 4. — Spectra of injected protons (dotted curves), protons surviving when the blob reaches 
the edge of the BLR (thin sohd curves), neutrons that escape from the blob (thick solid curves), 
neutrons (dashed curves) and gamma-rays (dot-dashed curves) that reach the edge of the BLR 

without attenuation for the same overall parameters as used in Figure 1 for 5 = 10, but calculated 
assuming different sizes of the blob, corresponding to t^ar = 0.3 and 3 days in the upper and 
lower panels, respectively. In calculation here we assume also that the same amount of accelerated 
protons as in Figure 1 is now injected continuously during the entire time of passage of the jet 
through BLR, 

tinj — 7.9 d in the observer frame, instead of tinj — 2d in Figure 1. 
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Fig. 5. — Radiation flux produced in and escaping from the blob (full dots) following the electromag- 
netic cascade initiated by energetic electrons and gamma rays produced in photopion interactions 
for the case 6 = 10 and the model parameters as in Figure 1. The thick and thin curves correspond 
to synchrotron and Compton scattered radiation, respectively. The radiation of the first gener- 
ation of electrons, which includes both the electrons from tt^ decay and the electrons produced 
by absorption of 7r°-decay gamma rays in the blob, are shown by the solid curves. The dashed, 
dot-dashed and 3-dot-dashed curves show contributions from the 2d, 3d and 4th generations of 
cascade electrons, respectively. The upper and lower panels (a) and (b) show the cascade radiation 
flux at t = 2d, when injection of accelerated protons stops, and at t = 7.9 d, when the blob reaches 
the edge of the BLR, respectively. 
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Fig. 6. — Total energy spectra of injected protons (dotted curve), and protons remaining in the blob 
(thin solid curve) at the time when the jet, with model parameters relevant to the April 1997 flare 
from Mrk 501 (see Section 3.2), tyar = 1 d and S = 10, is at a distance R = 0.4 pc from the center. 
The thick solid curve shows the total spectrum of the escaping neutrons, and the dot-dashed curve 
shows the spectrum of gamma rays that could be produced by these neutrons outside the blob and 
escape the scattering region. 
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Fig. 7. — Energy spectra of protons neutrons, and gamma rays from a jet in MrkSOl type object, 
as in Figure 6, but calculated for the case t^ar = 0.1 d and S = 25, with B = B^q = 1.0 in (a), and 
B = 0.2Beq in (b). 
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Fig. 8. — Mean-free-path of gamma rays due to pair production interactions with photons of the 
diffuse cosmic microwave background radiation for 2; = (soUd curve) and z = 0.538 (dashed 
curve), including the diffuse extragalactic infrared radiation field (see Section 4.3). 



